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Abstract 

We present a simple two-dimensional model of a px + ipy superfluid in which the mass flow that 
gives rise to the intrinsic angular momentum is easily calculated by numerical diagonalization of 
the Bogoliubov-de Gennes operator. We find that, at zero temperature and for constant director 
1, the mass flow closely follows the Ish.ika,Wci-]V[errn.in-]VIuzikcir formula jmass — 

icurl(pM/2). 
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I. INTRODUCTION 



Controversies over the "instrinsic angular momentum" of the A phase of superfluid ^He 
have a long history , but the issue has contemporary relevance for the search for 

boundary currents in the suspected + ipy superconducting phase of Sr2Ru04 j^, and for 
experiments on possible p-wave versions of atomic Fermi condensates ^|. 

Behind the competing computations of the ground-state angular momentum density lie 
two competing intuitive pictures. In the first the superfluid is regarded as a Bose condensate 
of Cooper pairs, each pair possessing angular momentum h. In this picture, and at tempera- 
ture T = 0, the total angular momentum of a spatially uniform system of A/" particles should 
be ^Afh. The second picture begins with the observation that in the BCS regime, where 
A ep, the opening of the energy gap A affects only states lying within an energy range 
of a few A about the Fermi surface at E = ep. It is therefore anticipated that the "naive" 
estimate of ^Afh will be reduced by a factor of A/ep |5l], or even by a factor of (A/ep)'^ Q]. 

For a system consisting of a fixed (even) number of particles with a common, angular 
momentum = 1 pair wavefunction ?/'(ri, r2), the antisjTiimetric many-body wavefunction is 
given by the pfaffian Pi[ip{ri, Vj)], and is undoubtedly an eigenstate of the angular momentum 
operator with eigenvalue ^HAf The problem is that ip{ri, is usually computed for 
a spatially uniform system, and a uniform system is not suitable for computing angular 
momentum: a substantial contribution to the angular momentum can arise from a small 
current at large distances — in particular near the walls of the container, where the density 
abruptly drops to zero. 

Once we allow spatial variations of the fluid density or order parameter, analytic expres- 
sions can only be obtained by approximate methods. Unfortunately, different approaches 
have led to different answers (see ref. jl| for a brief summary). The Gorkov gradient expan- 
sion, for example, computes Green functions for the Bogoliubov-de Gennes (BdG) equation, 
and in this formalism the angular momentum arises from the spectral asymmetry of the 
BdG differential operator. This asymmetry is closely related to the phenomenon of topolog- 
ically induced fractional charge, and to the axial anomaly of QED. Now the axial anomaly 
is quite a subtle phenomenon and can easily be missed if one makes naive manipulations 
in conditionally convergent integrals, or appeals to symmetries that are vitiated by bound- 
ary conditions. We claim however, that when solved exactly, the BdG equation produces 
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results entirely consistent with the Cooper-pair wave function approach: there is no A/e^ 
suppression, and the ground-state intrinsic angular momentum is per particle. 

In this paper we use the BdG equation to obtain numerical results for the angular momen- 
tum and associated mass-flows in a two-dimensional model of px + ipy superfiuid fermions 
confined in a harmonic trap. Numerical computations of the angular momentum in a cylin- 
drical container have been obtained by Kita , and our results are entirely consistent with 
this earlier work. We believe, however, that the simplicity of the present model, where the 
phenomena can be studied interactively with a few lines of Mathematical^ code, make it of 
interest. We begin in sections HTl and UTTl with a brief review of the Bogoliubov de Gennes 
formalism and the two-dimensional Harmonic oscillator. In section IIVI we show how sym- 
metries and the harmonic oscillator selection rules serve to reduce the BdG equation to a 
tridiagonal matrix eigenvalue problem which is solved numerically in section |Vl We end in 
section with a brief discussion of what we conjecture to be the origin of the much smaller 
angular momentum found by some older methods. 

II. SUPERFLUIDITY AND THE BOGOLIUBOV-DE GENNES EQUATION 

Because of the history of controversy in this subject, it is essential that we explain exactly 
how we do our calculations. We begin, therefore, with a brief review of the Bogoliubov-de 
Gennes formalism. 

Suppose that Hij is an A^-by-A^ matrix representing a one-particle {i.e. first quantized) 
Hamiltonian H. The second-quantized many-body hamiltonian corresponding to this is 

H = dlHijaj, (1) 

where a sum over the repeated indices is to be understood. Here a| and Sj are second- 
quantization fermion creation and annihilation operators obeying 

{cLi, ttj} = {a], a]} = 0, {oj, a]} = 5ij. (2) 

The operator a| creates a particle in state i and destroys such a particle. If we represent 
if in a continuous basis, the index "z" should be understood to incorporate both the space 
co-ordinate x, and any spin index. A sum over i therefore implies both an integral over real 
space and a sum over spin components. 
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To account for the effect of a condensate of Cooper pairs, we allow pairs of particles to 
disappear into or appear out of the background, and the second-quantized H is replaced by 
the Bogoliubov Hamiltonian 

-^Bogoliubov = alHijdj + jAijaloj + |A|^ajaj. (3) 

The gap-function matrix Aij is skew symmetric. Different forms of Ajj describe different 
condensates and result in different patterns of symmetry breaking. The entries in Aij are 
usually determined by a self consistency condition, or gap equation. In the present work 
we are more concerned with the consequences of a non-zero Aij than with its origin, and so 
will take its magnitude and form as something imposed externally. 

We can write the particle-number non- conserving Bogoliubov Hamiltonian as 

i^Bogoiiubov = 2 ( " V AJ . -H^- ) ^ 

where denotes the transpose of the hermitian matrix H. The two-by-two block-matrix 
form of the many-body Hamiltonian conveniently allows it to be diagonalized by means of a 
Bogoliubov transformation. We first solve the single-particle Bogoliubov-de-Gennes (BdG) 
eigenvalue problem 

/ H A \ f u„.\ , , 

Em\ . (5) 

\Vm J 

Here Um and Vm are A^-dimensional column vectors, which we take to be normalized so that 
l^mp -|- = 1. If we explicitly display the column-vector index i, these vectors can be 
regarded as A^-by-A^ matrices with entries Uim and Vim- Taking the complex conjugate of 
(IS) tells us that 

(6) 

and so the BdG-operator eigenvalues come in ± pairs. We will always use the symbol 
to refer to the positive eigenvalue, and a sum over otherwise unspecified eigenvectors is 
a sum over the positive eigenvectors. The manner in which the BdG operator doubles 
the spectrum, and how we retain only the positive part, is illustrated in Figure [H 
We now set 

4 = ViJ}^ + u*J}^^. (7) 






FIG. 1: The BdG operator spectrum for a one-dimensional non-relativistic gas with 
E(k) = k'^ /2 — fi. Left: the spectrum with A = 0. Right: the spectrum after a gap is opened 
by a non-zero A. In each case the part of the spectrum with BdG eigenvalue ii^ > is shown as a 
solid line, and the part with E < is shown dashed. For A 7^ 0, the E > "particle-like" region 
has ~ and |m| ~ 1, and the "hole-like" region has |m| ~ and \v\ ~ 1. 

The mutual orthonormality and completeness of the eigenvectors {um, VmY' ensures that the 
hm-, hln have the same anti-commutation relations as the dj, dj. In terms of the 6^, iln-, the 
second-quantized Hamiltonian becomes 

N ^ N ^ N 

-^Bogoliubov = ^ Emh\j)m " " ^ Em + - ^ Ef\ (8) 

m=l m=l i=l 

Here the Ef are the eigenvalues of Hij. Unlike the these can be of either sign. If 
all the quasi-particle excitation energies Em are strictly positive, the new ground state is 
non degenerate and is the unique state |0)^ annihilated by all the hm- The ground-state 
expectation value of an operator O = a\Oijaj is therefore 

(6) = ,(0|alO,,a,|0), 

= b{^\{ViJ)ra + U*mPm)Oij{Ujnhn + V*J)i)\Q) 

= VimOijVjm 
N 

= J2^0^m- (9) 

m=l 

Because of the E ^ —E, u ^ v* symmetry, we could equivalently write this last expression 
as 'Yli'^mQum, the sum being taken over the N negative-energy eigenstates of -f^BogoUubov, 
which we think of as being a filled Dirac sea. 

For example, the average number of particles present in the system is found by taking 
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Oij = Sij, and is 

N 

M = b{0\aldi\0)^ = VimV*^ = J2 l^rnl"^ = Yl l^'"'^" ^^^^ 

m=l negative _E 

When A = 0, this sum will be equal to the number of negative-energy eigenstates of Hij. 
Although we will not make use of it, it is worth pointing out that if we define the matrix 

= '"inXu'^Tmj, (11) 

then then the orthogonality and completeness conditions tell us that ipij = —"^ji- The new 
ground-state |0)^ can be written in terms of tpij as 

|0), = Cexp{lala]^,,}|0),. (12) 

Here |0)^ is the vacuum state annihilated by all the cij, and C is a normalization constant. 
This result exhibits the ground state as a coherent superposition of paired states, with 
ipij being the un- normalized Cooper-pair wavefunction. The state |0)j, gives a non-zero 
expectation value for the fermion- number non-conserving order parameter 

b{0\aiaj\0)f^ = ^UinV*^, (13) 

n 

and so possesses a sharp order-parameter phase and an uncertain particle number. If we 
desire to work with a fixed (even) number of particles J\f, we should retain only the |A/'-th 
term in the series expansion of the exponential in equation (IT2l) . The resulting many-body 
wavefunction is then the pfaffian 

^(n, . . . , W) = a{0\K ■ = Pf[^...,]- (14) 

In the large A/" limit there should be no locally measurable physical distinction between the 
sharp-phase and the sharp-particle-number ground states. 

III. TWO-DIMENSIONAL HARMONIC OSCILLATOR 

We consider a model that consists of spinless fermions confined in a harmonic trap in 
the x-y plane. The one-particle Hamiltonian H of section [ITl is therefore that of the two- 
dimensional harmonic oscillator: 

H = l{pl + pl + uj^x' + y')). (15) 
6 




Here, for example, px = —ihd/dx and so [x,px] = ih. We have taken the mass of the trapped 
particles to be unity. We define the usual harmonic-oscillator ladder operators 

and similarly ay and a^. These operators obey 

[ax, 4.] = [ay, al] = 1, [a^, ay] = [aj., ay] = 0, (17) 

and, in terms of them, 

H = huj{al.ax + ayay + 1). (18) 
The normalized eigenstates can be written as 

\nx,ny)=^^^-4=\0,0), (19) 



and have energy eigenvalues En^^ny = hw^rix + Uy + 1). 

We will find it more useful to work in a basis in which the angular momentum operator 

Lz = xpy - ypx = ih{a\ax - a^ay) (20) 

is diagonal with eigenvalues hi. To construct this basis, we define new ladder operators 

A = + ^2 = ^(^^^ " ^4)' (21) 

which obey 

[Lz,h\] = hh\, [Lz,bl] = -hbl (22) 

Consequently b\ increases the angular momentum quantum number / by unity and de- 
creases it by unity. It terms of these new operators we have 

H = huj{b\bi + b\b2 + 1), (23) 

and the eigenstates become 

\n,l)= ^ l—^ ^ |0,0), n = N + M, I = N - M. (24) 

The angular momentum of the state \n, I) is M and its energy is = ^(n + 1). The set 
of states with energy quantum number n is (n + l)-fold degenerate, the angular-momentum 
quantum number / running from / = —n to Z = n in steps of two. (See Figure [2j) 
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FIG. 2: The low-lying part of the harmonic oscillator spectrum. The zig-zag paths indicate how 
states are coupled by the tridiagonal matrix H^''\ For / positive, the first entry in the eigenvector 
will be a , and the second entry a "n," and so on. For / negative, the n's and u's are interchanged. 

We will require the normalized real-space wavef unctions (r, 6\n, I) corresponding to these 
eigenstates. These wavefunctions are commonly written as 



^^,i{r,e) = {r,e\2N + = 



7r(Ar+ |/|)! 



;^'^r'e--^'/2^i;iK2), (25) 



where l\^\x) is the associated Legendre polynomial 



X) 



-\m\ X Jn 



(26) 



which is of degree n. As the notation (r, 9\2N + |/|, /) indicates, these states have energy 



(27) 



and so the integer N counts the height (starting from zero at the lowest point) of the state 
in the column of eigenstates of angular momentum /. 



IV. BOGOLIUBOV-DE GENNES EIGENSTATES 

We now suppose that some suitable interaction has caused our fermions to enter a super- 
fluid phase characterized by an order parameter with + ipy symmetry. We can think of the 
fluid as being a single atomic layer of ^He in the A-phase and with the angular momentum 
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BdG 



(28) 



ium,i{r,0) 









(30) 



director 1 of the Cooper pairs pointing in the +z direction. We therefore wish to diagonahze 
the Bogohubov-de Gennes operator 

H - fi A{p^ + ipy) 
A{p^-ipy) -{H - fi) 
Here H is the harmonic oscillator hamiltonian, /i is a chemical potential that controls the 
number of particles in the trap, and A is a scalar parameter. The off-diagonal term 

A(p^. + ipy) = iAVn^{b\ - 62) (29) 

increases the angular momentum quantum number of any state on which it acts by unity, 
and represents the effect of the condensate of Cooper pairs, each pair possessing angular 
momentum +hz. 

The eigenstates with eigenvalue E = E^^i will be of the form 

The sum over n is over all harmonic-oscillator states consistent with the / quantum number. 
The pair (ra, I) corresponds to the index "i" in the general theory in section [TTl We choose to 
label the BdG eigenstates by the angular momentum of the lower component "t;." Because 
of the unit off-set between the angular momentum of the m's and f 's, the E ^ —E pairing 
is between eigenstates states with angular- momentum label / and those with label —I — 1. 

In the harmonic-oscillator eigenstate basis, and for any fixed /, the BdG operator reduces 
to a tridiagonal matrix connecting only the I and / + 1 subspaces. The pattern of ifedG- 
connected states is shown in Figure[2]for both positive and negative /. This pattern depends 
on the sign of / because, for I is negative, the ladder operator 62 appearing in px + ipy acts 
non-trivially on ||/|,— to yield a state with lower energy. When / is positive, however, 
it annihilates The factor of i in the upper entry of the column vector has been 

inserted to eliminate Vs and minus signs in the tridiagonal matrix, so that it becomes a real 
symmetric matrix with positive off-diagonal terms. 

The necessary matrix elements can be read-off from 



i{h\-h2)\n,l) = 



n + l + 2 



n + l,l + 1) -i 



n 



I. 



\n 



!,/ + !) 



i{bl - bi)\n,l + 1) = i 



{H-fj,)\n,l) = en\n,l), 



n-l-1. ^ . n + l + 1. ^ 
\n + 1, 1) — i\ \n — 1,1} 



(31) 



where 



e„, = hw^n + 1) — /i. 

For / positive or zero, and with huj = 1, the {I, I + l)-subspace BdG matrix takes the form, 

■ -ei AVTTT 
AVT+I q+1 Av^ 



-e/+2 AVr+2 
A^/^T2 Q+3 AV2 
Ay/2 -ez+4 



(32) 



All entries more than one step away from the diagonal are zero. When / is strictly negative 
the (— — 1/| + l)-subspace BdG matrix becomes 

■ e_i_i Av^ 
A^ AVI 

AVT t-i+i Av/^ITT 

A^-l + 1 -e_;+2 Ay/2 

AV2 e.i+3 ■ • . 



(33) 



In either case, the i, j-th matrix element links states with i,j labeling the distance from the 
lowest point along the zig-zag paths shown in Figure 1. The entries in the m-th eigenvector 
of H^^^ will therefore consist of alternating u'^/s and v^^/s, the first entry being a "f" when 
/ is positive, and a "m" when / is negative. 

It is an easy task to numerically diagonalize a tridiagonal matrix, and so the effect of the 
gap parameter A on the spectrum can be explored. 



V. NUMERICAL RESULTS 

We have computed ^ the spectrum and the eigenvectors and for a variety of 
values of fi and A. Typical results are displayed in this section. In all these plots we have 
set hw to unity. Changing the value of u serves only to rescale the energy and r. 

Figure [3] shows a plot of the eigenvalues for the case /i = 40.1 and A = 0.5. The 
interpenetrating wedges of ±E copies of the harmonic oscillator spectrum are clearly visible. 
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FIG. 3: Part of the BdG operator eigenvalue spectrum for // = 40.1 and A = 0.5. We have 
suppressed the axes for clarity. The horizontal co-ordinate, the angular momentum eigenvalue I, 
runs from —55 to +55, and the eigenvalues E^.i are plotted vertically. Each column of eigenvalues 
is the result of diagonalizing a 100-by-lOO tridiagonal matrix. 

as is the gap lying symmetrically about about e = 0. The family of states crossing the gap 
from the upper "continuum" to the lower as / increases is a chiral Majorana edge mode 
whose existence was first pointed out by Volovik , and whose importance for the analogy 
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with the Pfaffian quantum Hall state was stressed by Read and Green 
We also show a plot (Figure Sj) of the fluid density 

pir) = J2Mr.0)\' (34) 

1,171 

as a function of radius r. The A = 0.5 density profile differs from that for A = only by 
being slightly smoother. It closely follows the Thomas- Fermi density 

PTF(r) = < (35) 
to, r > 

The Thomas-Fermi approximation estimates the particle number to be 

ATtf = 2tt JJ^ p^F{r)rdr = ^/i^. (36) 

11 



p (r) 



\ 
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FIG. 4: Fluid density p{r) as a function of radius for the case fj, = 40.1, A = 0.5. The number of 
particles in the fluid is M = 2-k p{r) rdr = 840. 
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FIG. 5: Angular momentum density L{r) for the same parameters as Figure HI The total angular 
momentum is Ltot = ^tt Jq°° L{r) rdr = 410. Thus, for these parameters, Ltot/AA = 0.49. 

The actual particle number for the free-particle case A = is given by Aq = | [/ij ( [/ij + 1), 
where [/iJ indicates the integer part of /i. When A becomes non-zero the particle-hole 
asymmetry in the harmonic oscillator density of states causes the particle number at fixed 
/i to creep upwards. The effect is small for small A, however. 
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FIG. 6: The ground state angular momentum per particle L^^^/N plotted as a function of A for 
11 = 40.1. The discontinuous steps in the rising part of the graph are due to individual levels 
crossing e = 0, and so changing their occupation. The A > 0.05 ground state is therefore not 
adiabatically connected to the A = ground state. 

Figure [5] plots the ground state angular momentum density 

L{r) = (4) = J2^m,i (-^^) = -J2l\^^Ar,0)\' (37) 

l,m l,m 

as a function of r for the same parameters, and Figure [6] shows the angular momentum per 
particle as a function of A. Except for very small A we find that 

1 



2iT I L{r) rdr ^ -M. 



(38) 



There is no simple identity lying behind this fact, and mathematically it results from a quite 
non trivial rearrangement of spectral weight between the positive and negative E eigenstates 
of any given /. Physically, however, the angular momentum arises from an azimuthal mass 
flow 
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masSjO 



2i ^ 

Lm 



ld_ 

r 86 



ld_ 

r do 



-L(r) 



(39) 



This quantity is plotted in Figure [71 The straight line mass-flow distribution is not fun- 
damental, but indicates that the mass current is well described by the Ishikawa-Mermin- 
Muzikar formula [l2, 3, 3 

1 . ^l 



jmass = ^Cm\ ( -ph\ 



(40) 
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FIG. 7: The ground state azimuthal mass flow Jmass,e(^) = L{r)/r corresponding to the angular 
momentum density in Figure [5j 

and that the fluid has a near-parabolic particle density profile. The Friedel-like oscillations 
near the abrupt drop of Jmass,6»('") to zero at the edge of the droplet of confined fluid are also 
seen in the numerical results of ref. [8|]. 

If the mass current jmass is indeed given by ( HOi) . then an integration by parts gives 

(41) 



Ltot = Jir X j,„ass) d^r = j pldh = ^Afhl 



The quantity ^pHl can therefore be identified with an intrinsic angular momentum density. 
The current associated with this angular momentum density is then analogous to the bound 
electric current Jbound = curl M associated with a magnetic-moment density M. 

It is interesting to ask how much of the mass flow and angular momentum is supplied 
by the chiral Majorana edge mode. This branch of states is, after all, the most strikingly 
I ^ —I asymmetric feature of the BdG spectrum. A flrst (but misleading) estimate suggests 
that these states account for the entirety of the angular momentum. Figure [3] shows that 
positive-energy within-gap states exist for each integer / in the range —fi to 0. Because they 
are confined near the fluid boundary by Andreev reflection, each of these states consists of 



an equally- weighted linear superposition of particle and hole, and so has \u\ 
They therefore contribute an angular momentum of 

^edge-n.ode = -i ^ = 1^(;, + l)H ^ ^Mh = L^of 
l=—ii 



1 

2' 
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FIG. 8: The coeffecient \v\'^ for the lowest positive energy modes. The point at which the edge- 
modes merge into the upper continuum is signalled by the sharp decrease in \v\'^ near I = —40. 
The parameters are the same as those in Figure HI 

This result should be contrasted with what happens in Kita's model of a uniform fluid 
bounded by a rigid wall [8j]. In this case the edge modes have dispersion E{k) = —A{k/kp) 
2l\ and merge into the continuum at k = kp. They therefore contribute a boundary 

current — equivalently a momentum per unit length — of magnitude 



■ edge— mode 
.^boundary 



1 



nkz 

2 2n 



h 

Svr 



kp 



h 



-Pbulk, 



where pbuik = kp/An is the bulk fluid density. This current is larger by a factor of two than 
the boundary current jboundary = ^Pbuik/4 obtained from fHOl) by making p decrease from Pbuik 
to zero as we approach the container wall. The rigid-wall edge modes therefore oversupply 



angular momentum by a factor of two. It was shown in [21| that this twice-too-large bound 
state angular momentum is reduced by contributions from the unbound continuum states, 
and that the resulting edge momentum density is exactly what is required to give the N'h/2 
total angular momentum. 

To investigate whether the edge modes are the source of the entire harmonic trap edge 
current, we have isolated their contribution to the angular momentum density, and to the 
mass- flow current. Figure [S] shows that it is easy to determine which of the low- lying positive 
energy states states should be considered Andreev bound-state edge modes. When we sum 
the contributions to the angular momentum density and the mass flow from these states 
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FIG. 9: The contribution of edge modes with —41 < / < to the angular momentum density L(r) 
The parameters are those of Figure HI 
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FIG. 10: The —41 < Z < edge mode contribution to the mass flow. The parameters are the same 
as those in Figure HI 



only, we obtain the results shown in Figures [9] and [lOl Although 

POO 

(r) rdr 2n L{r) rdr 



27r / ^edge-mode/ 
'0 



(42) 



as anticipated, both L'^'^^^ raode^^^ JmassfT'^'^'^i^) localized near the boundary of the 
fluid, and differ substantially from L(r) and jmass,6»('^)- We conclude that, as with the rigid 



wall model of |8|, 



2l|, the continuum modes provide an important component of the mass- 
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flow. That the bound-state angular momentum contribution turns out to be equal to the 
total angular momentum should therefore be regarded as a coincidence arising from the 
particular form of the harmonic trap density profile. 



VI. DISCUSSION 

There are a number of gradient-expansion results for the mass current in three- 



dimensional ^He-A. By mapping the problem onto one involving fractional charge 15|, and 



for a uniform mass mass-density po, Garg et al. 



IGj obtained 



jmass = PoVs + ^Po curl HI - ^po 1(1 " curl h\) (43) 

where is the superfiuid velocity. This expression coincides with that obtained under the 

same conditions in 6(]. Because it does not allow for variations in the density, it cannot be 

j 

used to compute the boundary currents. Mermin and Muzikar [L'^ used a more sophisticated 
gradient expansion and found that when p is allowed to vary slowly on the length scale of 
the coherence length we have 

jmass = pvs + ^curl pkl - ^cq 1(1 ■ curl HI). (44) 

Here cq is a number that in the BCS limit A ^ ei;' is close to p, but goes to zero in the 
limit of tightly-bound Cooper pairs. The last term in both (H3l) and (jH]) is known as the 
the twist term. The current associated with it is now understood to be history dependent: 
if we start in the ground state with I spatially constant and then adiabatically deform 1 to 
the desired texture, the twist-term contribution to the mass flow is cancelled by momentum 
carried by excitations that have been forced through the nodes of the gap by spectral flow 
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(a manifestation of the QED axial anomaly) 

For a spatially constant director field 1, and when Vg = , t 
reduces to the earlier formula obtained by Ishikawa et al. jisl . |lj] : 



m 1191, 12L, 

le Mermin-Muzikar result 



jmass = ^CUrl Qp^l^ • (45) 

As we described in section |Vl although it is derived only for slowly varying p, our numerical 
results fit the two-dimensional version of fHS]) rather well. In particular, the striking straight- 
line dependence of jmass,6»(^) is an immediate consequence of the inverted-parabola Thomas- 
Fermi density profile being a good fit to the actual density. 
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Since many computations of the boundary current and resulting angular momentum 
are in agreement in their common domain of applicability, why do other estimates of the 
angular momentum, such as 3, Bk find results that are suppressed by powers of A/ej?? One 
explanation is that one can cast [19 



2l| the problem of the computing the boundary current 



j{k,) =lim|-i^sgn(K,feJ|^n,fcJ 'I 



into a weighted sum over of the quantity 

(46) 

where En^k^ the energy eigenvalues of a Dirac hamiltonian 

-f^Dirac = -icr-sdx + 0^2/^^ + m(x)(Ti (47) 
in which m{x) changes sign as x passes through zero. Now the operator 

Q = CTa i^Dirac " h (48) 

obeys 

Q-f^Dirac = — -f^DiracQ, (49) 

and so it seems that all the eigenstates of //Dirac, except for the E = topologically bound 
state which is anihilated by Q, come in pairs, ip and Qip, with opposite energy. If this 
pairing were really valid, then all terms, again with the exception of the unpaired bound 
state, would cancel in fj46|) . The formal E ^ —E symmetry suggested by the existence 
of Q is illusory, however. In order for ifoirac to possess a well defined eigenvalue problem 
we must impose some self-adjoint boundary conditions on the eigenfunctions. If ip obeys 
these boundary conditions, then, in general, Qip will not, and the actual density of states is 
asymmetric across E = 0. If, as was the case in the days when were written, one does 
not know of the topologically bound edge state, or is unaware of the fatuous nature of the 
symmetry implied by operators such as Q, then it might seem that the only contributions 
to the edge current come from the 0{A/eF) particle- hole asymmetry that arises from the 
curvature of the dispersion relation near the Fermi surface — exactly as claimed in ^,0]. 

In conclusion, we see that when solved exactly the BdG formalism produces a mass flow 
and angular momentum that coincides that obtained from the Cooper-pair wave function: 
there is no A/ep suppression, and the ground-state intrinsic angular momentum is per 
particle. 
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